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This  document  contains  three  technical  papers: 

1.  "Surface  Flaw  Effects  on  Crack  Propagation  from  Subsonic  Water  Drop  Impact" 

The  effects  of  surface  flaws  on  ceramic  target  damage  from  subsonic  water  drop 
impacts  are  investigated  using  nuncrical  code  simulations  based  on  fundamental  wave 
propagation  and  fracture  mechanics  concepts,  vj'he  water  drop  impact  generates  a 
tensile  stress  wave  with  a steep  near- surface  gradient"!  As~a'”consequence  of^this 
steep  stress  gradient,  qrack  activation  is  dependent  on  the  depth  of  surface  flaws 
In  addition,  a relatively  large  flaw  can  perturb  the  stress  field  such  that  crack 
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20. v activation  in  its  vicinity  is  suppressed. 

2.  "Effects  of  Surface  Pores  in  Silicon  Nitride  Impacted  by  a Water  Drop" 

A water  drop  impacting  directly  on  a surface  pore  may  significantly  increase 

the  probability  of  crack  formation  and  growth  near  the  pore  boundary . _ 

examine  this  physical  process,  the  dynamics  of  a 335  m/s  water  drop  impact  on 
a silicon  nitride  target  with  a 40  micron  diameter  hemispherical  surface  pore 
has  been  simulated  using  a finite  difference  computer  code  called  WAVE.-I~_ 

The  maximum  tensile  stress  in  the  silicon  nitride  target  is  about  15  kbar  = 

1.5  GPa  and  occurs  at  the  bottom  of  the  hemispherical  pore  on  the  axis  of 
cylindrical  symmetry.  The  tensile  stresses  at  the  bottom  of  the  pore  have  a 
duration  of  about  15  ns. 

Microstructural  failure  of  the  silicon  nitride  material  near  the  bottom  of  the 
pore  appears  likely  based  on  the  large  predicted  tensile  stresses  and  the 
X^large  predicted  stress  intensity  factors. 

3.  "Strain  Energy  Release  Associated  with  Circumferential  Cracks  Around 

Embedded  Particles" 

Numerical  code  simulations  are  used  to  determine  the  strain  energy  release 
associated  with  circumferential  cracks  around  spherical  particles  embedded  in 
tension  within  homogeneous  matrix  materials.  ,^It  is  shown  that  the  fractional 
release  of  strain  energy  f (p) , where  u is  the  normalized  crack  size  and 
f(0)=l  for  no  crack  and  f Ctt) =0  for  a crack  all  around  the  particle,  has  a 
point  of  inflection  at  about  p=tt/2  (hemispherical  crack) . A single  point  of 
inflection  in  f(y)  confirms  the  theoretical  conditions  of  crack  extension  and 
arrest  proposed  by  F.  F.  Lange  in  "Criteria  for  Crack  Extension  and  Arrest 
in  Residual,  Localized  Stress  Fields  Associated  with  Second  Phase  Particles" 
(Fracture  Mechanics  of  Ceramics,  Vol.  2,  ed.  by  Bradt,  Uasselman  and  Lange, 
pp.  599-613,  Plenum  Press,  New  York,  1974). 
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SUMMARY 


This  report  contains  a collection  of  technical  papers 
which  cover  three  important  areas  in  the  development  of 
physical  relationships  between  microstructure  and  fracture 
in  ceramics.  Each  section  in  this  report  covers  one  of 
the  following  investigations: 

1.  Surface  flaw  effects  on  crack  propagation  from  subsonic 
water  drop  impact. 

2.  Effects  of  surface  pores  in  silicon  nitride  impacted  by 
a water  drop. 

3.  Strain  energy  release  associated  with  circumferential 
cracks  around  embedded  particles. 


SURFACE  FLAW  EFFECTS  ON  CRACK  PROPAGATION  FROM  SUKSONIC  WATER  DROP  IMPACT1 


Y.  Marvin  Ito,  Marlin  Rosenblatt 


California  Research  6 Technology , Inc 
Woodland  Hills,  Ca.  11367 


The  effects  of  surface  flaws  on  ceramic  target  damage  from  subsonic  water 
drop  impacts  are  investigated  using  numerical  code  simulations  based  on  funda- 
mental wave  propagation  and  fracture  mechanics  concepts.  The  water  drop  impact 
generates  a tensile  stress  wave  with  a steep  near-stir  face  gradient.  As  a con- 
sequence of  this  steep  stress  gradient,  crack  activation  is  dependent  on  the 
depth  of  surface  flaws.  In  addition,  a relatively  large  flaw  can  perturb  the 
stress  field  such  that  crack  activation  in  its  vicinity  is  suppressed. 
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1.  Introduction 

The  damage  of  a ceramic  target  from  a subsonic  water  drop  impact  is  dependent  on  the 
. tensile  stresses  generated  by  the  impact  and  the  criteria  for  the  extension  of  flaws  ex- 
posed to  the  tensile  stress  field.  Since  the  target  response  is  essentially  elastic,  the 
maximum  tensile  stresses  are  the  radial  stress  components.  Also,  the  tensile  stress  field 
. has  a steep  near-surface  gradient.  Thus,  these  radial  stresses  tend  to  produce  circumfer- 
ential cracks  which  initiate  at  the  surface  and  propagate  into  the  target  [l] 2 . 

As  a consequence  of  the  relatively  steep  near-surface  stress  gradient,  it  is  expected 
[ lj  that  crack  activation  is  dependent  on  the  depth  of  surface  flaws  (exposed  to  this 
stress  gradient).  It  is  the  intent  of  this  paper  to  analyze  this  effect  for  specified  im- 
pact conditions. 


2.  Approach 

A series  of  numerical  simulations  are  performed  in  order  to  examine  the  effects  of 
surface  flaws  on  zinc  selcnide  window  damage  from  subsonic  water  drop  impact.  The  normal 
impact  conditions  of  a 2KU  diameter  water  drop  onto  the  window  target  at  velocity  VG  is 
shown  in  figure  1.  A finite-difference  computer  code,  WAVE-L,  with  an  analytic  pressure 
loading  model  is  used  to  predict  the  target  response. 

WAVF.-I.  is  a two-dimensional  (axial  or  nlano  symmetry),  Lagrangian,  explicit,  finite- 
difference  code  based  on  the  HEMP  scheme  [ 2]  which  integrates  the  governing  partial  dif- 
ferential equations  of  motion  for  arbitrary  dynamic  problems  in  solid  and  fluid  mechanics. 
It  has  been  applied  In  a number  of  particle  impact  investigations.  Including  both  sub- 
sonic [jl]  and  hypersonic  [4j  erosion. 

The  analytic  pressure  loading  model  has  been  developed  for  simulating  subsonic  water 
drop  impacts  and  is.  discussed  elsewhere  |f>].  The  impact  model  specifies  the  loading  pres- 
sure on  the  target  surface  as  a function  of  radius  from  the  impact  point  and  time,  P(r,t). 

‘Research  supported  by  Office  of  Naval  Research  under  Contract  No.  N00014-77-C-0790. 

Figures  in  brackets  indicate  the  literature  references  at  the  end  of  this  paper. 


Tlie  dove  1 opment  of  this  analytic  loading  model  was  based  on  detailed  WAVE-L  numerical  cal- 
culations [5]  of  spherical  water  drop  impact  onto  a rigid  surface.  The  rigid  surface 
approximation  is  valid  when  the  target  surface  deflection  is  negligible  compared  to  the 
drop  radius.  For  elastic  materials,  this  condition  is  satisfied  if  the  target  material 
impedance  (density  times  wave  speed)  is  much  greater  than  the  water  impedance.  The  zinc 
selenide  targets  described  in  this  paper  have  an  impedance  which  is  14  times  greater  than 
water.  Hence,  the  decoupled  approach  is  justified  in  the  present  case. 

The  loading  pressure  function,  l’(r,t),  for  a water  drop  of  2-mm  diameter  (2R0)  impact- 
ing at  222  m/ s (730  fps)  velocity  (Vn)  onto  a zinc  selenide  (7-nSe)  window  target  has  been 
determined  ^5].  Figure  2 gives  the  maximum  impact  pressure  profile.  Due  to  the  spherical 
geometry  and  water  material  properties,  an  off-axis  pea k pressure  occurs  (at  early  times) 
which  is  about  twice  the  Hugoniot  value.  This  peak  pressure  occurs  when  the  velocity  of 
the  radius  of  contact  (r  in  fig.  1)  approximately  equals  the  speed  of  sound  (C0)  in 
water  [_5]  . 


This  impact  pressure  l’(r,t)  is  used  as  the  loading  boundary  condition  for  the  axisym- 
metric  WAVE-L  numerical  calculation  ot  the  ZnSe  target  response.  Figure  3 shows  the 
initial  computational  grid  for  the  target  material  using  80  cells  per  water  drop  diameter 
(2R0  = 2mm).  In  the  WAVL-l.  code,  particle  velocity  is  defined  at  the  corners  of  the  com- 
putational cells  while  tin-  stress  state  is  associated  with  the  center  of  the  cells. 

Figure  4 indicates  typical  elect:  ■ reeperute  of  a ZnSe  window  in  terms  of  the  particle 
velocity  field  and  principal  in-plane  (>•,::>  feticide  stress  field  at  0.1  ps.  The  velocity 
vectors  show  the  magnitude  and  direction  of  the  velocity  of  the  particle  mass  located  at 
the  tail  of  the  vector.  A lr>-m/s  (SO- fps)  scale  bar  is  shown  in  the  upper  right  corner 
ot  the  figure.  The  principal  tensile  stresses  are  indicated  by  linos  which  show  their 
magnitude  and  direction  in  the  plane.  (Hoop  tensile  stresses  are  not  indicated  on 
these  plots.)  A 1-bar  (14500-psi)  scale  bar  is  shown  in  the  upper  left  corner  of  these 
plots.  Thus,  at  0.1  ps  after  the  impact,  the  peak  velocities  are  roughly  15  m/s  (50  fps) 
in  the  target  and  the  peak  tensile  stresses  are  about  2 khar  (30000  psi). 

It  is  seen  that  a region  of  high  tensile  stresses  occurs  near  the  surface,  outside 
the  contact  area.  The  material  directly  under  the  contact  area  is  in  pure  compression. 
Figure  5 shows  the  corresponding  velocity  and  in-plane  tensile  stress  fields  at  0.2  ps. 

The  qualitative  features  are  similar  to  the  0.1-ps  plot;  however,  the  peak  velocities  and 
peak  tensile  stresses  have  already  decayed  below  the  0.1-ps  values. 

The  specified  impact  conditions  used  in  tlie  present  numerical  simulations  involve 
222  m/s  (730  fps)  impact  of  2-mm  diameter  water  drop  onto  zinc  selenide  windows.  These 
conditions  are  considered  to  be  representative  of  subsonic  liquid  drop  impact  onto  ceramic 
targets.  In  the  following,  Section  3 gives  details  of  the  near-surface  stress  gradient, 
Stction  4 shows  the  effect  of  surface  flaws  and  Section  5 has  the  concluding  discussion. 


3.  Near-Surface  Stress  Gradient 

The  formation  of  high  tensile  stresses  near  the  surface  of  a ZnSe  target  impacted  by 
a water  drop  is  illustrated  in  figures  4 and  5 (Section  2).  In  this  numerical  simulation 
the  initial  computational  grid  (fig.  3 in  Section  2)  has  40  cells  per  water  drop  radius. 
Since  the  effective  prediction  of  crack  activation  requires  that  these  high  tensile 
stresses  be  well  characterized,  an  appropriate  fine-zoned  computational  grid  must  be 
selected.  Some  calculations,  employing  selected  fine-zoned  computational  grids  in  the 
near-surface  impact  region,  provide  .1  basis  for  selecting  the  grid  structure  for  subse- 
quent crack  activation  experiments. 


'Note  that  this  pressure  peak  will  not  occur  for  subsonic  colid  particle  impacts; 
and  observed  differences  in  target  response  from  solid  versus  liquid  projectile 
impacts  may  be  partially  due  to  the  presence  of  the  off-axis  pressure  peak. 


ilO 


Figure  2.  Maximum  impact  pressure  as  a function  of  radial  distance  from  impact  center  for 
222-m/s  impact  of  2-rnrn  water  drop  onto  ZnSe  window. 
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Figure  5.  Elastic  response  of  ZnSe  window  for  222-m/s  impact  of  2-mm  water  drop 
at  0.2  us  (80  computational  cells  per  drop  diameter). 


Figure  6 shows  the  peak  tensile  stress  versus  depth  in  a ZnSe  window  impacted  by  a 
water  drop  at  222  m/s  (7 JO  fps)  for  three  variations  in  computational  cell  size.  Note 
the  very  steep  stress  gradient  near  the  impact  surface  and  its  increased  resolution1*  with 
decreased  cell  size  (where  stress  is  defined  at  the  cell  center).  For  a water  drop  of 
1-mm  radius5*  the  finest-zoned  (80  cells  per  drop  radius)  case  is  adequate  to  resolve  peak 
tensile  stress  from  356  MPa  (52000  psi)  at  a depth  of  about  6 pm  to  80  MPa  (12000  psi)  at 
100  pm  depth. 

Figure  7 shows  the  peak  tensile  (radial)  stress  near  the  surface  versus  radius.  The 
results  are  from  the  surface  computational  cells  which  have  stress  defined  at  the  cell 
center  depth  of  z * R,,/l60  or  6.25  pm.  The  maximum  in  the  peak  tension  (356  MPa)  occurs 
at  radius  rp  a 0.2  R0  or  200  pm.  For  r < 0. 1 R0  or  100  pm  no  tension  develops  in  the 
target  surface.  At  radius  2 rp  ~ 0.4  K„  or  400  pm  the  peak  tension  has  dropped  to  only 
200  MPa  (30000  psi).  The  nominal  static  tensile  strength  of  ZnSe  is  about  45  MPa 
(6500  psi). 

The  temporal  development  of  the  radial  tensile  stress  at  four  radial  locations 
(3/4  rp,  rp,  3/2  rp  and  2 rp)  near  the  target  surface  is  given  in  figure  8.  The  pulse 
duration  (which  is  about  30  ns  at  the  maximum  tension  location  rp)  increases  with  radial 
location  ('60  ns  at  3/2  rp  and  -75  ns  at  2 rp) . 

The  effect  of  a surface  flaw  of  depth  Ro/80  or  12.5  pm  at  radius  rp  = 0.2  R0  or 
200  pm  on  the  peak  ra< Hal  tensile  stress  is  shown  on  figure  9.  For  this  flaw,  the  stress 
gradient  ahead  of  the  crack  tip  is  similar  to  that  from  the  surface  without  a flaw,  but 
the  maximum  in  the  peak  radial  tension  is  about  10  percent  less. 

It  is  concluded  that  the  effects  of  surface  flaws  of  depths  from  10  pm  and  greater, 
which  arc  exposed  to  these  stress  gradients,  should  be  well  characterized  in  subsequent 
••rack  activation  experiments  using  the  fine-zoned  comput at ionu 1 grid  with  80  cells  per 
drop  radius. 


4.  Sui  t ace  Flaw  Kffocts 

In  the  numerical  code  a flaw  is  simulated  as  "tensile"  cracks  across  computational 
cells  near  the  surface.  In  a cell  with  an  initial  crack,  no  tensile  stress  across  the 
crack  is  permitted,  and  no  shear  stress  is  permitted  on  the  crack  if  the  crack  is  open 
(the  width  of  the  crack  is  continuously  monitored).  Those  stress  adjustments  are  impor- 
tant because  a crack  frequently  alters  the  local  stress  field  in  such  a way  as  to  enhance 
its  own  growth.  Thus,  any  realistic  method  of  predicting  crack  growth  must  consider  the 
altered  stress  field. 

The  crack  activation  experiments  involve  surface  flaws  of  various  depths  located  at 
radius  3/2  rp  a 0.3  Ro  or  300  pm.  At  this  location,  without  a flaw,  the  maximum  radial 
tensile  stress  is  about  285  MI’a  and  a tensile  stress  gradient  develops  to  a depth  of 
/.  ■ 0.075  R0  or  75  pm  (see  fig.  10).  Also,  the  tensile  pulse  duration  near  the  surface  is 
about  60  ns  (see  fig.  8 in  Section  3). 

Figure  10  shows  the  effect  of  surface  hoop  (circumferential)  flaws  on  the  radial 
tensile  stress  gradient  ahead  of  the  crack  tip  for  crack  depths  of  25,  50,  and  100  pm 
(H0=lmm).  In  each  case,  the  peak  radial  tensile  stress  ahead  (one-half  a cell  dimension) 
af  the  crack  tip  is  enhanced,  but  decreases  with  increasing  flaw  depth. 


'*  In  fact,  an  analytic  solution  by  Blowers  [ b]  predicts  infinite  peak  tensile  stress 
at  the  impact  surface  for  purely  elastic  targets. 

’ Note  that  these  elastic  solutions  can  |>e  scaled  linearly  In  dimensions  and  time  for 
other  drop  radii  as  long  as  the  impact  velocity  and  material  properties  are  not 
changed . 


ZnSe  window 


■/t- 


Figure  11  shows  the  effect  of  the  flows  at  radius  '1/2  Tp  ~ 0.3  R0  or  300  pm  on  the 
peak  radial  tensile  stress  gradient  monitored  downstream  at  radius  2 Tp  - 0.4  R0  or  400  pm. 
A factor  of  four  reduction  (from  the  no  flaw  case)  Is  observed  near  the  surface  for  the 
100  pm  (R0=lmm)  deep  flaw.  (Note  that  the  stress  profile  upstream  at  Tp  ~ 0.2  Rn  or 
200  pm  remains  essentially  unaffected  hv  the  flaws.) 


5.  Discussion 

The  effects  of  surface  flaws  on  ceramic  target  damage  from  subsonic  water  drop  imp. lets 
are  investigated  using  numerical  code  simulations.  The  pressure  pulse  due  to  water  drop 
impact  generates  a tensile  stress  wave  with  a steep  near-surface  gradient.  Figure  12 
shows  a schematic  view  of  the  effect  of  a surface  flaw  on  the  near-surface  tensile  stress 
profile.  Upstream  profiles  remain  essential lv  unchanged,  while  downstream  profiles  have 
reduced  gradients  and  peak  value  at  surface. 

Using  the  stress  intensity  approach  of  fracture  mechanics,  the  stress  intensity 
factor,  Kj,  can  be  obtained1',  as  shown  in  figure  12,  and  crack  activation  is  defined  by 
the  condition  where  K[  >.  Figure  13  shows  the  peak  stress  intensity  factor  for  various 

surface  flaws  based  on  the  results  of  the  crack  ju  t i vat  ion  experiments  from  Section  4 (see 
fig.  10).  For  a toughness  value  ot  Ep  = 1 Ml’a/m  , surface  flaws  of  depth  greater  than 
T)  pm  will  not  propagate. 

Thus,  as  a consequence  of  the  steep  stress  gradient,  surface  crack  activation  under 
water  drop  impact  is  dependent  on  the  depth  of  surface  flaws.  In  addition,  a relatively 
large  surface  flaw  can  perturb  the  stress  field  such  that  crack  activation  in  its  vicinity 
is  suppressed. 


The  authors  wish  to  acknowledge  the  useful 
University  of  California,  Berkeley  (formerly  at 
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^ The  stress  intensity  factor  is  calculated  from 


Kt  - o/2*!  = o/itA 

where  0 is  the  radial  stress  located  one-half  a cell  dimension,  A/2,  from  the  crack  tip. 
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ABSTRACT 


A water  drop  impacting  directly  on  a surface  pore 
may  significantly  increase  the  probability  of  crack  form- 
ation and  growth  near  the  pore  boundary.  To  examine  this 
physical  process,  the  dynamics  of  a 335  m/s  water  drop 
impact  on  a silicon  nitride  target  with  a 40  micron 
diameter  hemispherical  surface  pore  has  been  simulated 
using  a finite  difference  computer  code  called  WAVE-L. 

The  numerical  simulation  solves  the  finite  difference 
analogs  to  the  physical  conservation  equations  of  mass, 
momenta,  and  energy  along  v.ith  the  constitutive  equations 
for  the  water  and  silicon  nitride. 

Shortly  after  impact,  stress  wave  dynamics  accelerate 
the  free  surface  of  the  water  near  the  axis  as  this  water 
flows  into  the  pore.  The  water  on  axis  impacts  the  bottom 
of  the  pore  with  a velocity  which  is  about  a factor  of  7 
higher  than  the  initial  impact  velocity.  This  increase  in 
the  water  flow  velocity  and  the  associated  lateral  flow  of 
water  into  the  pore  lead  to  a relatively  large  pressure 
pulse  in  the  pore  when  the  water  impacts  the  bottom  of  the 
pore.  A peak  compressive  pressure  of  ~60  kbar  = 6 GPa  is 
predicted;  this  pressure  is  about  a factor  of  9 higher  than 
the  1-D  llugoniot  impact  pressure. 

The  maximum  tensile  stress  in  the  silicon  nitride 
target  is  about  15  kbar  = 1.5  GPa  and  occurs  at  the  bottom 
of  the  hemispherical  pore  on  the  axis  of  cylindrical  symmetry. 
The  tensile  stresses  at  the  bottom  of  the  pore  have  a duration 
of  about  15  ns. 

Microst rue tura 1 failure  of  the  silicon  nitride  material 
near  the  bottom  of  the  pore  appears  likely  based  on  the  large 
predicted  tensile  stresses  and  the  large  predicted  stress 
intensity  factors. 


EFFECTS  OF  SURFACE  PORES  IN  A S II,  I CON  NITRIDE 
TARGET  IMPACTED  BY  A WATER  DROP 


1.  INTRODUCTION 

A water  drop  impacting  directly  on  a surface  pore 
may  significantly  increase  the  probability  of  crack 
formation  and  growth  near  the  pore  boundary.  To  examine 
this  physical  process,  the  dynamics  of  a 335  m/s  water 
drop  impact  on  a silicon  nitride  target  with  a 40  micron 
diameter  hemispherical  pore  has  been  calculated.  This 
pore  size  is  representative  for  reaction  bonded  Si3N4 
(Reference  1).  The  pore  diameter  is  small  compared  to 
millimeter  size  water  drops  of  interest;  thus,  the 
water  boundary  is  assumed  to  be  a flat  surface,  as  indi- 
cated in  Figure  1.  The  stress  wave  dynamics  between  the 
water  drop  and  silicon  nitride  target  arc  simulated  using 
a finite  difference  computer  code  called  WAVE-L.  WAVE-L 
is  a Lagrangian  code  of  the  HEMP  type  (Ref.  2).  The 
numerical  simulation  solves  the  finite  difference  analogs 
to  the  physical  conservation  equations  of  mass,  momenta, 
and  energy  along  with  the  constitutive  equations  for  the 
water  and  porous  silicon  nitride.  The  water  is  treated  as 
a compressible  inviscid  fluid  while  the  silicon  nitride  is 
treated  as  linear  elastic.  The  material  properties  for 
these  materials  are  indicated  on  Figure  1;  the  detailed 
equation  of  state  for  water  can  be  found  in  Reference  3. 

The  basic  physical  response  of  the  water  and  target 
material  near  the  pore  is  described  in  Section  2.  The  dyna- 
mic tensile  stress  characteristics  in  the  target,  which 
provide  the  driving  forces  for  crack  initiation  and  growth, 
are  described  in  Section  3.  In  Section  4,  the  conclusions 
are  presented  and  discussed. 
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2.  THE  WATER  DROP  AND  TARGET  DYNAMICS  NEAR  THE 
SURFACE  PORE 

Figure  2 shows  the  initial  computational  grid  config- 
uration. The  left  axis  is  an  axis  of  cylindrical  symmetry. 

The  explicitly  defined  hemispherical  surface  flaw  is  treated 
as  a void  in  this  numerical  simulation.  The  computational 
cells  have  been  zoned  to  provide  a good  definition  of  stresses 
near  the  boundary  of  the  surface  flaw.  Stresses  are  calcu- 
lated at  all  cell  centers,  while  particle  velocities  are 
calculated  at  all  cell  corners.  Note  that  the  water  can 
slide  over  the  surface  of  the  porous  silicon  nitride  target. 
The  sliding  interfaces  are  considered  frictionless.  The 
sliding  occurs  both  along  the  hemispherical  boundary  of  the 
surface  pore  as  well  as  on  the  flat  front  surface  of  the 
target  near  the  pore. 

Figure  5 shows  a particle  velocity  (left  half)  and 
principal  stress  (right  half)  field  plot  at  8 nanoseconds  (ns) 
after  impact.  The  velocity  vectors  and  three  principal 
stresses  are  plotted  for  both  the  water  and  Si3N4  directly 
from  the  numerical  simulation.  Scale  liars  for  the  magnitude 
of  the  vectors  and  stress  components  are  indicated  at  the  top 
of  the  field  plots.  Also,  a description  of  the  principal 
stress  conventions  is  shown. 

The  velocity  flow  field  at  8 ns  shows  the  downward  and 
lateral  flow  of  water  into  the  pore.  The  lateral  flow  into 
the  pore  is  a direct  consequence  of  the  large  stress  gradients 
which  form  in  the  water  upon  impact.  Upon  impact,  the  1-D 
llugoniot  stresses  (P  = 6.8  Khar  = 680  MPa)  develop  near  the 
interface  between  the  water  and  silicon  nitride  target. 

However  in  the  region  of  the  pore,  the  water  does  not  impact 
any  target  material  and  therefore  no  pressures  initially 
develop.  The  initial  pressure  gradient  at  r = 20  pm  in  the 
water  causes  pressure  waves  to  propagate  toward  the  axis  of 
cylindrical  symmetry  and  causes  inward  radial  velocities  to 
develop  in  the  water. 


The  very  high  wave  velocities  in  the  silicon 
nitride  target  cause  the  entire  pore  boundary  to  be  stressed 
very  early  (-3  ns).  Note  that  by  8 ns  stress  waves  in  the 
target  have  propagated  to  depths  of  about  .008  cm  = 80  pm 
which  is  relatively  large  compared  to  the  pore  depth  of 
20  pm. 

Figure  4 shows  the  velocity  and  principal  stress  field 
at  15  ns.  Note  the  relatively  large  volume  of  water  which 
has  now  achieved  inward  radial  velocities.  The  peak  free 
surface  water  velocity  in  the  pore  is  now  82'  higher  than 
the  initial  impact  velocity.  In  the  silicon  nitride  target, 
the  particle  velocities  arc  very  low  due  to  the  much  higher 
impedance  (pCo)  in  the  target  as  compared  to  the  water. 

Figure  5 shows  the  corresponding  plot  at  23  ns.  By  this  time, 
the  water  along  the  boundary  of  the  pore  has  moved  over  half 
way  down  the  pore.  The  peak  particle  velocity  of  the  water 
is  now  ISOlhighcr  than  the  initial  impact  velocity. 

Figure  b shows  the  velocity  and  principal  stress  field 
plots  at  30  ns  after  impact  with  the  target  surface.  At 
this  time,  the  water  has  almost  completely  filled  t lie  pore 
with  only  a small  volume  near  the  pore  bottom  remaining  as 
void.  (Notice  that  the  effects  of  compression  of  the  air 
initially  filling  the  pore  have  been  ignored  in  this  calcula- 
tion). The  velocities  near  the  axis  are  now  quite  large,  a 
factor  of  -4  greater  than  the  initial  velocity.  Note  that 
both  the  velocity  and  stress  scale  bars  have  been  changed  in 
this  plot  due  to  the  large  spatial  field  of  view.  The 
stresses  in  the  target  near  the  pore  arc  in  tension  at  the 
bottom  of  the  pore  at  this  time.  In  the  next  section  these 
tensions  will  be  shown  to  reach  -IS  kb  = 1500  MPa  * 1.5  GPa. 

At  33  ns  (Figure  7)  the  water  impacts  the  bottom  of  the 
pore  and  a strong  compressive  pulse  is  generated  in  the  water 
and  in  <he  target  near  the  bottom  of  the  pore.  The  impact 


velocity  of  the  water  on  the  pore  bottom  is  a factor  of  ~7 
greater  than  the  initial  impact  velocity.  The  peak  compres- 
sive pressure  is  ~60  kbar  = 6, DUO  MPa.  The  pressures  propagate 
through  the  water  in  the  pore  and  begin  to  reverse  the  flow  of 
the  water  as  shown  on  Figure  8,  t = 38  ns.  By  53  ns  (Figure  9), 
water  is  flowing  out  of  the  pore. 

3.  MAXIMUM  PRINCIPAL  TENSILE  STRF.SS  HISTORIES 
IN  THE  POROUS  SILICON  NITRIDE  TARGET 


During  the  numerical  simulation,  several  target  locations 
were  monitored  for  principal  stress  histories.  Figures  10 
through  16  show  the  maximum  principal  tensile  stress  histories 
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Figure  10  shows  the  tension  history  at  Station  1 at  the 
bottom  of  the  pore  on  the  axis  of  cylindrical  symmetry.  (A 
sketch  appears  on  each  figure  showing  the  location  of  the 
monitored  stresses.)  The  kev  reference  time  is  t = .033  ps 
* 33  ns  which  is  when  the  water  impacts  the  bottom  of  the 
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pore.  The  greatest  tension  calculated  anywhere  in  the  target 
is  15  kbar  = 1.5  GPa  which  develops  at  Station  1 shortly  after 
water  impacts  the  bottom  of  the  pore.  However,  the  tensile 
stress  history  at  the  bottom  of  the  pore  shows  large  tensile 
stresses  of  nearly  the  same  magnitude  forming  just  prior  to 
t = 33  ns.  Then,  when  the  water  impacts  the  bottom  of  the 
pore,  a compressive  stress  wave  begins  to  propagate  into  the 
target  and  the  principal  stresses  near  the  bottom  of  the  pore 
are  all  compressive*.  Once  this  compressive  stress  pulse 
propagates  and  attenuates  in  the  target,  the  relatively  high 
pressures  present  in  the  water  filling  the  pore  cause  large 
"hoop"  tensile  stresses  near  the  bottom  of  the  pore.  These 
tensile  stresses  have  a duration  of  about  15  ns. 

Figure  11  shows  the  maximum  tensile  stress  history  at 
Station  2 (y  = 26°  near  the  pore  surface).  Also,  on  this 
figure  and  the  other  tensile  stress  histories,  the  Station  1 
results  are  shown  for  reference.  The  tensile  stresses  do  not 
get  as  high  or  last  as  long  off  the  axis  of  cylindrical 
symmetry  as  on  the  axis. 

Figure  12  shows  the  tensile  stress  history  at  Station  3, 

T = 49°.  Figure  13  shows  the  tensile  stress  history  at 
Station  4 , = 90°.  No  significant  tension  develops  at  this 

station . 

Figures  14,  15  and  lb  show  the  tensile  stress  histories 
at  axial  Stations  7,  9 and  10,  respectively.  note  that  the 
peak  tension  decays  rapidly  with  depth.  By  Z = 2.4U,  the  peak 
tension  is  reduced  to  ~2  kbars  = 200  MPa  and  has  a duration  of 
~ 2 ns . 

Figure  17  shows  a comparison  of  the  numerical  simulation 
stress  predictions  on  the  pore  surface  at  t = 45  ns  with  a 
simple  "steady  state"  analytical  model.  This  comparison  was 


* The  value  plotted  in  Figure  10  is  the  least  compressive 
principal  stress. 
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performed  to  provide  partial  validation  of  the  numerical 
simulation  results  for  this  configuration.  The  analytical 
model  is  composed  of  a superposition  of  two  solutions  found 
in  Timoshenko  (Ref.  4,  pages  392-398).  The  two  solutions 
are  combined  to  attain  a uniaxial  strain  far  field  boundary 
condition.  The  uniform  internal  pressure  component  was 
selected  to  approximate  the  water  pressure  in  the  pore. 

The  comparisons  are  very  good,  especially  near  the  axis  of 
symmetry  (4*  = 0)  where  the  relatively  high  water  pressure 
in  the  pore  dominates  the  formation  of  the  hoop  tension  (o^) 
and  radial  compression  (o  ) . 

4.  CONCLUSIONS 

Based  on  the  numerical  simulation  results,  the  following 
conclusions  arc  presented  for  a hemispherical  surface  pore  in 
silicon  nitride  impacted  at  335  m/s  by  a water  drop: 

1.  The  maximum  tensile  stress  in  the  silicon  nitride 
target  is  about  15  kbar  = 1.5  GPa  and  occurs  at  the 
bottom  of  the  hemispherical  pore  on  the  axis  of 
cylindrical  symmetry.  The  tensile  stresses  at  the 
bottom  of  the  pore  have  a duration  of  about  15  ns. 

2.  Stress  wave  dynamics  accelerate  the  water  near  the 
axis  as  this  water  flows  into  the  pore.  The  water 
on  axis  impacts  the  bottom  of  the  pore  with  a velo- 
city which  is  about  a factor  of  7 higher  than  the 
initial  impact  velocity.  This  increase  in  the  water 
flow  velocity  and  the  associated  lateral  flow  of 
water  into  the  pore  lead  to  a relatively  large 
pressure  pulse  in  the  pore  when  the  water  impacts 
the  bottom  of  the  pore.  A peak  compressive  pressure 
of  ~60  kbar  = 6 GPa  is  predicted;  this  pressure  is 
about  a factor  of  9 higher  than  the  1-J  llugoniot 
impact  pressure. 

3.  Microstructural  failure  of  the  silicon  nitride  material 
near  the  bottom  of  the  pore  appears  likely  based  on  the 
following  two  considerations.  First,  there  is  a re- 
ported static  tensile  strength  of  79-141  MPa  (Ref.  1) 
which  is  an  order  of  magnitude  lower  than  the  predicted 
dynamic  tensile  stresses.  Second,  the  "effective"  stress 
intensity  fadtor  exceeds  the  measured  Kjc.  Note  that  a 

. stress  intensity  factor  cannot  be  rigorously  defined  for 
the  ideal  hemispherical  pore,  since  there  is  no  crack  tip. 
However,  the  following  sketch  and  first  order  relation- 
ship represent  an  approximation  for  the  effective  stress 
intensity  factor  on  a small  crack  located  at  the  bottom 
of  the  pore. 
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I f we  assume  a small 
crack  at  the  bottom  of 
the  pore,  then  the  calcu- 
lated tensile  stresses, 
o,  a distance  -p  below  / | 
the  pore  bottom  are  re-  / | 
lated  to  the  stress 
intensity  factor,  Kj,  by  I* 
K i - o/lv'Zp 11)  Zp 


Using  o = 1500  MPa  and  Z^  = 2 pm  in  Equation  1 to 

obtain  an  approximate  value  for  the  maximum  predicted 

stress  intensity  factor,  we  find  K™ax  - 5 MPa-m’"  which 

is  greater  than  the  reported  critical  stress  intensity 

1 

factor  of  Kjp  = 2-4  MPa -in'*  for  reaction  bonded  S i 3 N «, 
(Ref . 1). 


The  largest  calcu- 
lated tensile  stress 
of  1500  MPa  occurred 
at  a depth  of  Zp  = 2 pm 


A numerical  simulation  approach  to  crack  initiation, 
propagation  and  arrest  for  water  drop  impacts  on 
ceramics  can  be  found  in  References  5 and  6.  This 
general  dynamic  crack  growth  approach  should  be 
applied  to  predict  crack  growth  characteristics  near 
a surface  pore  subjected  to  water  drop  impacts. 
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figure  2. 


Computational  Crid  Configuration  at  Initial  Condition. 
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Figure  4.  Particle  Velocity  and  Principal  Stress  Fields  at  1 5 Nanoseconds 
for  an  1100  fps  (335  m/s)  Water  Prop  Impact  on  a 40pm  Diameter 
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40um  Diameter  Pore  (1Kb  = lOOMPn). 
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SECTION  III 

Strain  Energy  Release  Associated  with  Circumferential 
Cracks  around  Embedded  Particles 


STRAIN  HNEIRfiY  RFLHASF  ASSOCI  ATl-I)  WITH  C I RCUMFF.RENTI  AL 
CRACKS  AROUND  LMBLDDF.l)  P ART  I CI.HS  1 


Y.  Marvin  lto,  Martin  Rosenblatt,  Louis  Y.  Cheng 

California  Research  l\  Technology,  Inc. 
Woodland  Hills,  Ca . 91367 


Numerical  code  simulations  are  used  to  determine 
the  strain  energy  release  associated  with  circumfer- 
ential cracks  around  spherical  particles  embedded  in 
tension  within  homogeneous  matrix  materials.  It  is 
shown  that  the  fractional  release  of  strain  energy  f(u), 
where  u is  the  normalized  crack  size  and  f (01=1  for  no 
crack  and  f(n)=0  for  a crack  all  around  the  particle, 
has  a point  of  inflection  at  about  u=n/2  (hemispherical 
crack).  A single  point  of  inflection  in  f(u)  confirms 
the  theoretical  conditions  of  crack  extension  and 
arrest  proposed  by  Lange  1 l]2. 


I n t roduc t i on 

An  important  area  for  developing  new  ceramic  materials 
involves  second  phase  particles  embedded  within  a matrix  phase. 
Due  to  the  usual  difference  in  thermal  expansion,  highly  local- 
ized stresses  develop  as  the  ceramic  is  cooled  from  its  fabrica- 
tion temperature.  Thus,  these  inclusions  are  usually  considered 
(1)  as  possible  precursors  to  cracks  that  form  either  during 
cooling  or  during  subsequent  stressing.  The  mechanisms  which 
govern  the  formation  of  these  cracks  arc  therefore  a critical 
part  of  the  micro- st ructural  behavior  of  these  ceramics. 

For  the  case  of  a single  spherical  inclusion  of  radj_us  R 
within  an  infinite  medium,  a uniform  hydrostatic  stress  o 
develops  within  the  particle  and  radial  and  circumferential 
stresses  or=o  R J/ r 3 and  op=-o  R3/2r3  develop  within  the  surround 
ing  matrix  rapidly  falling  off  with  increasing  distance  from  the 
inclusion  as  shown  on  Figure  1.  The  stress  o due  to  temperature 
change  AT  (positive  for  cooling)  is  given  (2]  by 

— AaAT  , , 


1 Research  supported  by  Office  ol  Naval  Research  under  Contract  N00014-77-C-0790. 

Figures  in  brackets  Indicate  the  literature  references  at  the  end  of  this  papei 
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Normalized  Stress  in  Matrix 
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where  k ~ (l  + vm)/2Hm  ♦ l 1 - 2\>pi /lip  ami  A - 1 <tp-um,  I-  is  Young's 
modulus,  v is  Poisson's  ratio,  and  u is  thermal  expansion  co- 
efficient. Also,  tension  is  positive  and  subscripts  m and  p 
denote  matrix  and  particle,  respectively.  (It  should  be  noted 

that  the  stress  o is  independent  of  particle  size.)  The 
associated  strain  energy  is  [ 3) 


2 11  Sj’iJ-AVpJ  ]<3 

■p 

, ( 1 vm)  ..j 

“7;  k 

m 


- 2n  k o*  R 3 (2) 


Now  consider  a crack,  which  is  specified  by  a single 
dimensional  parameter  c,  associated  with  the  spherical 
particle  and  favorably  oriented1  to  t lie  high  tensile  stress 
field  near  the  particle/matrix  interface.  The  strain  energy 
associated  with  the  crack/particle/matrix  interaction  for  a 
given  normalized  crack  size  u = c/R  can  be  expressed  (1)  as 

USE  = usi: 

where  f(p)  is  a dimensionless  function  which  defines  the 
fractional  strain  energy  release  and  by  definition, 

l^f (y)>o. 

By  hypothesizing  that  the  strain  energy  release  f(p) 
possesses  an  inflection  point,  Lange  ( 11  has  determined  the 
general  criteria  for  crack  extension  and  arrest  in  the  local- 
ized stress  fields  associated  with  second  phase  particles. 

The  unique  result  of  this  analysis  is  that  the  particle  size  is 
one  of  the  key  factors  that  governs  the  criteria.  Thus,  for  a 
given  material,  a critical  size  exists  below  which  crack  exten- 
sion cannot  occur  and  is  given  by 


R 


cons  t ant 

o ■’ 


(4) 


3 


For  particles  In  tension  (Aa^O)  cracks  can  occur  either  within  the  particle 
or  within  the  matrix.  For  particles  in  compression  (Aa<0)  cracks  should 
only  occur  within  the  matrix. 


This  is 
13]  . 


diiect  agreement  with  expo  r i mont  a I observations 


It  is  tlio  intent  ol  this  paper  to  determine  the 
explicit  function  for  the  release  of  strain  energy  associ- 
ated with  circumferential  cracks  within  the  highly  localized 
stress  field  due  to  a spherical  particle  embedded  in  tendon 
within  a homogeneous  matrix  material. 


Approach 

A c i rcum ( e r en t i a 1 crack  associated  with  a spherical  particle 
or  inclusion  placed  in  tension  within  a homogeneous  material  has 
the  following  problem  description: 


Axis  of 


S >1111110  tpy-xj 

J:  Mat 


Crack 


Crack. 


lH  r 


Angle 


U * — 

M R 


4 tir2 ( l-cos9) 


(Tin*  inclusion  Is  in  hydrostatic 
tension  and  the  matrix  Is  In  radial 
tension  and  hoop  compression.) 


A series  oi  nun, idea l simulations  are  performed  in  order  to 
examine  the  quantitative  effects  of  these  cracks.  The  simulated 
conditions  of  a single  spherical  particle  o I radius  R within  an 
infinite  matrix  is  shown  in  figure  2.  An  inner  core  region  of 
radius  ;|  = R/*1  and  an  outer  annulus  region  of  radius  b*3R  arc 
simulated  by  appropriate  pressure  boundary  conditions1*.  A finite- 
difference  computer  code,  WAVP-I,,  with  the  analytic  representation 
of  the  initial  stress  field  is  used  t o predict  the  variation  of 
the  stress  field  with  crack  size. 


The  strain  enoigy  associated  with  those  regions  is  less  than  5 percent 
of  the  total. 


■ 


WAY'li-1.  is  a t wo-d  imens  iona  1 (axial  or  plane  symmetry), 
Lagrangian,  explicit,  finite-difference  code  based  on  the  IIHMP 
scheme  [4]  which  integrates  the  governing  partial  differential 
equations  of  motion  for  arbitrary  dynamic  or  quasi -static  problems 
in  solid  and  fluid  mechanics.  It  has  been  applied  in  a number  of 
fracture  mechanics  investigations,  including  subsonic  water  drop 
impact  [5]  and  quasi -static  indentation  [ 6 ) of  ceramics. 

The  analytic  representation  of  the  highly  localized  stress 
field  (sec  Fig.  1)  is  used  as  the  initial  conditions  for  the 
axisymmctric  WAVF.-I.  numerical  calculations  of  the  crack/inclusion/ 
matrix  interactions.  Figure  3 shows  the  initial  computational 
grid  using  a cell  size  of  Ar  0.03R  near  the  inclusion/matrix 
interface  to  resolve  the  high  initial  stress  gradient  to  within 
10  percent.  In  the  WAVF.-l.  code,  velocity  (and  displacement)  is 
defined  at  the  corners  of  the  computational  cells  while  stress 
(and  strain)  state  is  associated  with  the  center  of  the  cells. 


In  the  numerical  code,  cracks  arc  simulated  as  "tensile 
failure"  across  computational  cells.  In  a cell  with  a crack,  no 
tensile  stress  across  the  crack  is  permitted,  and  no  shear  stress 
is  permitted  on  the  crack  if  the  crack  is  open  (the  width  of  the 
crack  is  continuously  monitored).  The  quasi  static  stress  (and 
strain)  field  is  progressively  calculated  as  a crack  is  increment- 
ally  extended  along  a prescribed  path  (at  the  particle/matrix 
interface  in  the  present  investigation). 


St  t a i n F.ne  rgy  Re  1 ease 

In  the  present  numerical  simulations,  the  strain  energy  re- 
lease is  calculated,  as  a circumferential  crack  is  extended 
along  the  particle/matrix  i n tor  face 5,  for  embedded  particles  placed 
in  tension,  for  three  modular  ratios.  Table  1 gives  the  material 
properties  and  initial  conditions  for  the  three  cases  with  modu- 
lar ratios  of  M= 1 , M 1/4  and  M l. 

Figure  I shows  the  calculated  fractional  strain  energy 
release  function  f(|i).  Due  to  accumulated  numerical  error  as  the 
crack  extends,  all  the  strain  energy  is  not  released  (f(n)40)  as 
the  crack  extends  all  around  the  particle  (u=n) . However,  the 
error  is  no  greater  than  10  percent  for  Case  2 (M=l/4)  and  5 
percent  for  Case  1 (MM)  and  Case  3 ( M - 4 ) . 

Within  this  numerical  accuracy,  the  present  calculations  show 
no  significant  effect  of  modular  ratio  on  strain  energy  release. 
In  all  three  cases  there  is  a single  point  of  inflection  in  f(p) 
at  about  u = ti/2.  This  is  shown  more  clearly  on  Figure  5 where  the 
derivative  df/dp  has  a minimum  value  at  about  p-tt/2  and  compares 
well  with  4 curve  fit  with  symmetry  about  p*it/2.  Note  that 
p=n/2  corresponds  to  a hemispherical  crack  and  may  be  primarily 
due  to  geometry  effects. 


^ The  crack  is  associated  with  [lie  computational  cells  in  the  matrix  next  to 
the  interface. 
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Figure  3.  Initial  Computational  Grid  for  Simulated  Problem  of  a Spherical  Particle  of  Radius 
Embedded  Within  Infinite  Matrix. 


TABLE  1.  Matt-rial  Properties  and  Initial  Conditions  for  Three  Modular  Ratios 


r Boundary  at 
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Figure  4.  Fractional  Strain  Energy  Release,  f(u).  Versus  Normalized 
Crack  Size,  V1  < for  Three  Cases  with  Modular  Ratios  of 
M * 1 , M ■ . and  M * 4 . 
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Figure  5.  Fractional  Strain  Energy  Derivative,  df/dp,  Versus 

Normalized  Crack  Size,  |i , and  Comparison  with  Analytic 
«r  Curve  Fit. 


An  analytic  curve  fit  to  the  calculated  results  is  given 
in  the  form 


( n _ df 
(u)  = 


- a s 1 n • h + 


sin  h 


and 


(5) 


f(u) 


— , U sinw  c o s vi  -> 
'a(T 2 1 


♦H( 


COS  3lJ 

3 


cosu)  + c 


lb) 


Using  the  conditions  f ( 0 ) h 1 , f(u)sO  and  f'(n/2)=0.55  (see 
Fig.  5),  the  three  constants  are 


a = 1.12 
E = 0.57 
c = 1 . 58 


(7) 


Figures  5 and  6 show  this  analytic  function. 


D i s c u s s i o n 


The  associated  strain  energy  release  rate''  is  defined  hv 

dUj.  j. 

G H - -as- 


where  A = 4 tt R 2 ( 1 -cosw)  is  the  surface  area  of  the  circumfer- 
ential crack.  It  follows  from  Equations  (2),  (3)  and  (5)  that 
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- U 


0 f ' ( H ) 

S H 4 ti  R 7 s i n vi 
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(9) 


Hence,  from  Lange  [1],  using  the  thermodynamic  criterion  intro- 
duced by  Griffith  (7],  crack  extension  will  occur  when  the  value 
of  G for  a given  crack  size  is  equal  to  y and  crack  arrest  occurs 
at  a value  of  \i  where  G-Y  , where  y is  the  fracture  energy.  This 
is  shown  in  Figure  7 using  Case  1 properties  (sec  Table  1)  and  a 
nominal  value  of  Y = 10  J/m2 . Regardless  of  the  crack  size,  crack 


* The  stress  intensity  factor  can  he  derived  through  K « /C*F.  '. 
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Figure  6.  Analytic  Representation  of  Fractional  Strain  Energy  Release, 
f(p),  Versus  Normalized  Crack  Size,  p. 
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Figure  7.  Strain  Energy  Release  Rate  for  Crack  Extension  and  Arrest 
Using  Case  1 Properties  from  Table  1. 


The  present  investigation  obtained  (i  for  the  case  of 
circumferential  cracks  associated  with  isolated  embedded 
particles  placed  in  tension.  Radial  cracks  [1]  are  observed 
for  embedded  particles  placed  in  compression,  and  thus  would 
require  a separate  solution  for  G.  Also,  future  solutions 
would  involve  part i cl e- t o-part ic 1 e interactions  for  the  case 
where  the  inter-particle  separation  is  less  than  one  particle 
diameter . 
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